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ABSTRACT  39<0       °L 


In  modern  turbofan  engines,  variable  geometry  has  been 
incorporated  to  improve  some  off-design  performance.  Most 
control  designs  ignore  this  variable  geometry  and  use  fuel 
metering  as  the  primary  control  input. 

This  thesis  investigates  the  use  of  variable  geometry  to 
control  the  engine  and,  thereby,  reduce  fuel  consumption  due 
to  transients.   Additionally,  steady-state  trim  conditions 
are  altered  to  reduce  the  static  fuel  consumption.   The  non- 
linear transient  simulation  program  is  used  to  analyze  the 
steady-state  operating  condition  and  develop  small  perturba- 
tion control  limitations.   Linear  models,  both  large  and 
reduced  order,  are  used  in  analyzing  the  effect  of  controllers 
on  system  response.   A  computer  program  was  generated  to  re- 
duce a  large  order  linear  model  to  a  usable  size  for  control 
system  development. 

This  analysis  shows  the  reduced-order  regime  dependent 
controllers  to  be  viable  and  to  favorably  enhance  the  quest 
for  reducing  specific  fuel  consumption  in  existing  engines. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION 11 

II.  TURBOFAN  ENGINE  PERFORMANCE  CONSIDERATIONS  14 

A.  FLIGHT  REGIME  FACTORS  I4 

B.  TURBOFAN  ENGINE  PERFORMANCE  MEASUREMENTS  17 

C.  ENGINE  CONTROL  DEVICES  19 

D.  ENGINE  MODELING 21 

E.  CONTROL  DEVICE  EFFECTIVENESS  25 

III.  CONTROL  DEVELOPMENT  34 

A.  LINEAR  OPTIMAL  CONTROL  THEORY  34 

B.  CONTROL  SYSTEM  MODELS  36 

C.  REDUCTION  OF  LINEAR  MODELS 4i 

D.  REDUCED  ORDER  MODELS  45 

E.  VERIFICATION  OF  REDUCED  ORDER  MODELS  52 

IV.  ANALYSIS  OF  CONTROL  SYSTEMS  55 

A.  COST  FUNCTION  WEIGHING  FACTORS  55 

B.  REDUCED  ORDER  SIMULATION  ANALYSIS  66 

C.  FULL-ORDER  SIMULATION  ANALYSIS  68 

V.  CONCLUSIONS  AND  RECOMMENDATIONS  73 

APPENDIX  A 75 

APPENDIX  B 95 

LIST  OF  REFERENCES  103 

INITIAL  DISTRIBUTION  LIST  105 


LIST  OF  TABLES 

II-l  Engine  States  24 

II-2  Engine  Inputs  and  Outputs 26 

II-3  Input  Perturbations  31 

II-4  Summary  of  Derivatives  ^2 

III-l  Operating  Point  Eigenvalues ^6 

III-2  Reduced  Model  Comparisons  4^ 

III-3  Reduced  Model  Eigenvalues  49 

III-4  Fifth  Order  Model  50 

III-5  Third  Order  Model  51 

IV-1  Case  Definition 58 

IV- 2  Actuator  Limitations  6  0 


LIST  OF  FIGURES 

II-l  Flight  Envelope  Limitations  15 

II-2  Flight  Regimes 16 

II-3  Engine  Cutaway 18 

II-4  Compressor  Stage  Flow 22 

II-5  Change  in  SFC  vs  Control  Inputs 28 

II-6  Change  in  Thrust  vs  Control  Inputs  29 

II-7  Change  in  Turbine  Inlet  Temperature  vs 

Control  Inputs  30 

III-l  Model  Comparison  in  Frequency  Domain  53 

IV-1  SFC  vs  Control  Inputs 62 

IV-2  Thrust  vs  Control  Inputs 63 

IV-3  Fan  Speed  vs  Control  Inputs 64 

IV-4  Compressor  Speed  vs  Control  Inputs  65 

IV-5  Third  Order  vs  Fifth  Order  Control  Response  67 

IV-6  Fifth  Order  vs  Full  Order  Control  Response  69 

IV-7  Third  Order  vs  Full  Order  Control  Response  70 

IV-8  Full  Order  Implementation 71 


TABLE  OF  SYMBOLS 

SFC  Thrust  specific  fuel  consumption 

M3  Total  engine  airflow 
a 

Mf  Main  burner  fuel  flow 

U  '  Hot  stream  velocity 
e  * 

U  •  •  Cold  (fan)  stream  velocity 

pn  Total  pressure  at  station  x 

Tn  Total  temperature  at  station  x 

R  Universal  gas  constant 

Fn  Net  engine  thrust 

A.  Nozzle  jet  area 

CIW  Inlet  variable  guide  vanes 

RCW  High  compressor  variable  vanes 

X  The  state  variable  vector 

Y  The  output  vector 

U  The  input  vector 

x$  The  $  element  of  the  X  vector 

y$  The  $  element  of  the  Y  vector 

u$  The  $  element  of  the  U  vector 

F  The  plant  matrix 

Fr  The  reduced-order  plant  matrix 

G  The  input  gain  matrix 

Gr  The  reduced-order  input  gain  matrix 

H  The  output  matrix 

Hr  The  reduced-order  output  matrix 


D  The  feed-forward  matrix 

Dr  The  reduced-order  feed-forward  matrix 

K  The  feedback  gain  matrix 

R,  The  state  weighting  matrix 

R~  The  input  weighting  matrix 

P  The  Riccati  solution  vector 

T  The  transpose  matrix 

Z  The  redefined  state  vector 

A  The  output  control  matrix 

A  The  output  weighting  matrix 

J  The  cost  function 

8  The  fan  bypass  ratio 

n  The  efficiency  factor 

Y  The  ratio  of  specific  heats 


ACKNOWLE  DGMENTS 

I  would  like  to  thank  all  of  the  people,  without  the 
help  of  whom,  this  thesis  could  not  have  been  completed. 

I  would  especially  like  to  thank  Professor  D.  J.  Collins, 
my  thesis  advisor,  for  his  help  and  guidance;  Mr.  Neil  Harvey 
of  the  W.  R.  Church  Computer  Center;  and,  most  especially, 
my  wife  and  daughter  for  their  support  and  patience  during 
my  tenure  at  the  Naval  Postgraduate  School. 


10 


I.   INTRODUCTION 

High  performance,  military  aircraft  are  operationally 
tasked  to  conduct  missions  throughout  their  designed  per- 
formance envelopes.   To  accomplish  these  missions,  the  air- 
craft must  operate  efficiently  within  its  flight  regime. 
The  typical  jet  engine  is  designed  for  operation  at  a  single 
specified  altitude  and  Mach  number  combination,  e.g.,  30,000 
feet  at  M  =  0.9.   This  design  point  will  give  inherently 
good  cruise  performance  and  is  particularly  well-suited  to 
commercial  transport  operations.   This  engine  design  method 
forces  the  user  into  accepting  off-design  performance  in 
typical  military  applications.   Engine  performance  needs  to 
be  optimized  to  efficiently  conduct  all  assigned  missions, 
e.g.,  air  combat  maneuvering,  ground  attack,  supersonic  inter- 
ception and  carrier  operations.   Adaptive  digital  control 
methods  can  be  incorporated  within  the  engine  to  accomplish 
this  mission  dependent  optimization. 

Modern  control  methods  allow  the  designer  to  develop  a 
control  system  to  regulate  every  parameter  of  the  parent 
system.  He  has  been  aided  by  computer  simulations  of  non- 
linear systems  and  linear  approximations  of  those  systems. 
The  desired  control  system  can  be  developed  parallel  to  or 
ahead  of  the  controlled  system. 

The  state-of-the-art  turbofan  engine  utilizes  variable 
geometry,  in  such  forms  as  variable  incidence  inlet  guide 
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vanes,  adjustable  exhaust  nozzle  area  and  airflow  bleed-off. 
Hy drome chan ical  controls  have  been  used  in  the  past  to  im- 
prove some  elements  of  engine  performance  through  positioning 
these  variable  engine  components.   Electronic  controls,  such 
as  were  introduced  in  the  Pratt  &  Whitney  F100-PW-100  turbo- 
fan  engine,  allow  improved  scheduling  of  the  variable  devices, 
again  improving  the  engine  performance.   Studies  have  been 
done  on  the  F100  engine  to  optimize  its  performance  at  a 
point  other  than  its  design  point  [1] . 

Originally,  the  onboard  microprocessing  capability  of  an 
aircraft  was  intended  for  weapon  system  improvement  and 
integration.   Further  computing  equipment  was  added  as  part 
of  the  engine  control  to  supervise  the  hydromechanical  con- 
trol system  and  detect  faults.   The  F100  engine  introduced  an 
electronic  control  which  allowed  use  of  digital  control  tech- 
niques incorporating  the  hydromechanical  devices  as  implemen- 
tors  rather  than  controllers.   Digital  device  improvements 
continue  to  expand  the  available  computing  power  available, 
making  flight  control  and  engine  control  integration  a  dis- 
tinct possibility.   As  more  software  space  becomes  available 
to  the  engine  controller,  regime  adaptive  control  will  be  the 
accepted  standard  in  military  applications. 

This  thesis  looks  at  using  modern,  multivariate  control 
logic  in  designing  flight-regime  dependent  control  for  a 
typical  military  turbofan  engine,  the  Pratt  &  Whitney  F100. 
Specifically,  the  attempt  is  made  to  improve  thrust  specific 
fuel  consumption,  i.e.,  fuel  economy,  in  the  high  altitude 
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cruise  environment.   Such  improvement  in  fuel  economy  is 
especially  important  in  the  conservation-conscious  decades 
ahead. 
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II.   TURBOFAN  ENGINE  PERFORMANCE  CONSIDERATIONS 

A.   FLIGHT  REGIME  FACTORS 

The  ability  of  a  military  aircraft  to  perform  an  assigned 
mission  is  determined  by  its  flight  regime.   Figure  II-l  is 
a  typical  flight  regime  for  a  modern,  high-performance 
military  aircraft. 

The  confines  of  the  flight  regime  are  dictated  by  numer- 
ous aerodynamic  and  propulsive  factors.   A  well  designed 
aircraft  will  incorporate  an  engine  having  operating  boundar- 
ies in  excess  of  those  of  the  airframe.   This  design  technique 
allows  for  a  factor  of  safety  in  actual  operation  as  well  as 
compensate  for  installation  and  equipment  aging  losses. 

At  the  low  speed,  low  altitude  region  of  the  envelope, 
engine  performance  is  judged  in  the  ability  to  produce 
sufficient  thrust  for  take-off  and  landing  and  to  provide 
rapid  engine  acceleration  in  the  missed  approach  or  balked 
landing.   The  extreme  upper  region  of  the  flight  envelope 
requires  the  engine  to  produce  maximum  thrust  while  not  ex- 
ceeding turbine  temperature  limitations. 

Figure  II-2  is  a  mission  delineated  flight  regime  for  the 
supersonic  tactical  fighter  aircraft.   The  desired  mission 
performance  from  this  chart  must  be  interrelated  with  the 
aerodynamic  and  propulsion  limitations  of  Figure  II-l.   An 
adaptive  control  system  can  be  designed  to  schedule  optimum 
performance  of  the  airframe  and  engine  throughout  the  flight 
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Figure  II-l.   Flight  Envelope  Limitations 
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Figure    II-2.      Flight   Regimes 
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regime.   Northrup  Corporation  is  involved  with  a  U.S.  Air 
Force  project  to  implement  an  adaptive  flight  control  system 
in  an  aircraft  such  as  the  F-16.   The  same  type  of  approach 
can  be  applied  to  engine  control  and  integrated  into  an 
aircraft  control  system  to  optimize  the  total  package. 

B.   TURBOFAN  ENGINE  PERFORMANCE  MEASUREMENTS 

The  primary  factor  with  which  to  judge  engine  performance 

is  the  engine's  ability  to  convert  fuel  to  thrust  efficiently 

This  measurement,  thrust  specific  fuel  consumption,  SFC,  is 

a  function  of  the  airflow  through  the  engine,  M  ,  the  bypass 

a 

ratio,  b,  and  the  fan  and  core  exhaust  velocities,  U  ' '  and 

e 

U  ',  also  termed  cold  and  hot  stream  velocities, 
e 

M 
SFC   =   -±[(l+f)u  '  +  3U  '  '  -  (1+8)V]         (II-D 

1 1  ^  G  w 

where  Mf  is  the  fuel  flow,  f  is  the  fuel-to-air  ratio  and  V 
is  the  flight  velocity.   Figure  II-3  shows  the  F100-PW-100 
engine  in  cross  section  with  reference  labels. 

The  engine  designer  fixes  the  airflow  and  bypass  ratio  of 
the  engine.   The  exhaust  velocities  are  functions  of  the  ex- 
haust temperatures  and  pressures,  which  are  controllable 
through  fuel  metering  and  variable  geometry. 


2Yl  pa   Yt-IA,   1/2 


1 


e 


^l1  06        P06 


■y  Pi/        ^-f^- 

uo"     =      i~  R  Tns   nil  -  (~-)Y     '  Y]  }  (II-3) 

e  y-i  ua  p08 
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Figure    II-3.      Engine    Cutaway 


where  y  and  Yi  are  the  respective  ratios  of  specific  heats, 
Tnfi  and  Tna  the  respective  exhaust  total  temperatures,  n  is 
the  exhaust  nozzle  efficiency,  and  p  /p0fi  and  p  /p0g  the 
respective  ratios  of  free  stream  ambient  pressure  to  dis- 
charge total  pressure. 

The  standard  bypass  fan  will  produce  a  fixed  pressure 
ratio  and  exhaust  temperature,  producing  an  exhaust  velocity 
in  the  cold  stream  in  excess  of  the  flight  velocity.   This 
increases  the  engine's  thrust  without  altering  the  hot  stream 
through  the  engine  core  (ignoring  energy  lost  in  driving  the 
fan) .   The  advanced,  low-bypass  fans  that  reintroduce  the 
cold  flow  into  the  augmentor,  as  in  the  F100,  increase  the 
augmentor  efficiency  and  during  non-augmented  operation, 
allow  for  increasing  U  ' '  through  the  use  of  the  variable  ex- 
haust nozzle. 

Studies  show  that  Thrust  Specific  Fuel  Consumption,  SFC, 
can  be  decreased  through  increased  turbine  inlet  temperature 
and  compressor  pressure  ratio.  Both  of  these  variables  are 
designed  into  the  engine  at  the  engine  design  point.  Engine 
control  hardware  can  be  designed  to  maintain  the  turbine  inlet 
temperature,  FTIT,  and  pressure  ratio,  R  ,  at  the  optimum  or 
best  achievable  level  during  off-design  operations. 

C.   ENGINE  CONTROL  DEVICES 

The  mechanical  devices  available  for  engine  control  are 
the  hydromechanical  fuel  control,  the  variable  exhaust  nozzle 
and  variable  incidence  inlet  guide  vanes  or  stators.   The 


19 


first  device  has  the  effect  of  increasing  the  fuel-to-air 
ratio  and  temperatures  in  the  turbine  and  exhaust  sections  of 
the  engine.   The  combined  effect  of  increasing  f  and  U  '  is 
increased" thrust,  but  does  not  necessarily  improve  SFC  due  to 
the  corresponding  increase  in  the  fuel  flow.   The  variable 
area  exhaust  nozzle  has  been  used  for  the  last  two  decades  in 
most  high  performance  military  aircraft.   This  device  in- 
creases the  hot  section  exhaust  velocity,  U  ' ,  thus  increasing 
thrust  and  maintaining  near  ideal  pressure  in  the  exhaust 
section  of  the  engine. 

The  latest  mechanical  devices  are  the  variable  incidence 
blading  concepts.   The  first  use  of  these  devices  was  to 
achieve  near  design  mass  flow  rate  through  the  compressor 
during  the  engine  start  cycle.   These  variable  stators  would 
be  rotated  to  a  starting  configuration  until  the  engine  reached 
a  given  rotational  speed  and  then  be  driven  to  a  fixed  normal 
operating  position. 

The  full  advantage  of  variable  geometry  blading  comes  in 
the  optimizing  of  off -design  performance.   The  single  stage 
pressure  ratio  is  found  to  be 

Pm  U   2  C 

-^     =      1   +    n    5^ — [1    -   n^<tan  bo    +    tan    ai>l  (H-4) 

p01  RT01     Ur      2        l 

where  n  is  the  stage  efficiency,  U   is  the  rotational  velocity 
of  the  blading,  R  is  the  Universal  Gas  Constant  and  TQ,  the 
total  temperature  at  the  front  face  of  the  stage.   At  the 


20 


design  condition,  the  factor  (tan  b~  +  tan  a,)  is  minimized, 
thus  providing  the  best  pressure  ratio.   The  angle  b~  is  the 
relative  flow  direction  upon  leaving  the  rotor  portion  of  the 
machine;  a,  is  the  angle  at  which  the  flow  impinges  the  rotor 
and  is  a  function  of  the  stator  angle  (see  Fig.  II-4) .   Con- 
trolling the  stator  angle  increases  the  stage  pressure  ratio 
in  the  early  stages  and  increases  the  compressor  (or  fan) 
performance  over  a  wide  range  of  operations. 

In  the  F100-PW-100  engine,  four  mechanical  devices  are 
used,  fuel  control,  variable  exhaust  nozzle,  inlet  guide  vanes 
with  moveable  trailing  edge  and  moveable  stators  in  the  first 
three  compressor  stages.   The  exhaust  nozzle  maintains  nozzle 
area,  expansion  ratio  and  boattail  drag,  simultaneously,  near 
optimum.   The  inlet  guide  vanes  are  used  ahead  of  the  fan  to 
improve  inlet  distortion  tolerance,  improve  fan  efficiency 
and  enhance  engine  acceleration  performance.   The  variable 
stators  improve  starting  and  high  Mach  number  characteristics. 

D.   ENGINE  MODELING 

It  is  not  economically  practical  to  build  an  engine  and 
then  conduct  experimentation  to  determine  its  performance. 
Modern  technology,  in  both  digital  and  hybrid  computers, 
allows  the  design  engineer  to  build  a  simulation  model  of  the 
engine  and  predict  the  system  response  well  in  advance  of 
hardware  manufacture.   Several  generic  computer  simulations 
are  in  use,  these  allow  non-specific  engine  performance  analy- 
sis.  Additionally,  a  computer  model  is  generated  for  each 
new  development  engine  by  the  engine  manufacturer. 
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Figure  II-4.   Compressor  Stage  Flow 
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Szuch,  [2] ,  provides  a  detailed  description  of  current 
engine  modeling  techniques.   The  Pratt  &  Whitney  F100-PW-100 
engine  has  been  modeled  with  both  a  real-time  hybrid  computer 
[3],  and  a  digital  transient  simulation  [4].   Both  programs 
have  distinct  advantages  and  applications.   Although  the  digi- 
tal program  is  not  real-time,  it  provides  a  large  volume  of 
data  and  is  easily  accessed  and  maintained.   A  copy  of  the 
F100-PW-100  Transient  Engine  Simulation  Deck  (CCD  1103.3.0) 
was  obtained  from  the  YF100  Special  Project  Office  of  the 
Aeronautical  Systems  Division  of  the  U.S.  Air  Force  at  Wright- 
Patterson  Air  Force  Base,  Ohio.   This  program  was  modified 
to  be  compatible  with  the  Naval  Postgraduate  School  IBM 
370/3033  computer  system. 

Designing  a  controller  based  on  the  non-linear  digital 
simulation  alone  is  extremely  difficult,  if  not  impossible. 
Miller  and  Hackney  [5] ,  developed  high  order  linear  approxi- 
mations, at  steady  state  operating  points,  from  the  non- 
linear simulation  program.   The  linear  system  was  modeled  in 
the  classical  state  variable  form 

X   =   FX  +  GU  (II-5) 

Y   =   HX  +  DU  (II-6) 

The  16  system  states,  listed  in  Table  II-l,  were  chosen  to 
coincide  with  the  states  measured  in  the  non-linear  simulation, 
These  states  include  2  rotational  speeds,  3  internal  total 
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TABLE  II-l 
Engine  States 

xl  —  Fan  speed,  (Nl) ,  in  RPM 

x2  —  Compressor  speed,  (N2) ,  in  RPM 

x3  —  Compressor  discharge  pressure,  (Pt3),  in  psia 

x4  —  Interturbine  volume  pressure,  (Pt4.5),  in  psia 

x5  —  Augmentor  pressure,  (Pt7m) ,  in  psia 

x6   —   Fan  inside  diameter  discharge  temperature,  (Tt2.5h), 
in  Rankine 

x7   —  Duct  temperature,  (Tt2.5c),  in  R 

x8   -  Compressor  discharge  temperature,  (Tt3)  ,  in  R 

x9   —  Burner  exit  fast  response  temperature,  (Tt4hi) ,  in  R 

xlO  —  Burner  exit  slow  response  temperature,  (Tt41o) ,  in  R 

xll  —  Burnerexit  total  temperature,  (Tt4) ,  in  R 

xl2  —   Fan  turbine  inlet  fast  response  temperature,  (Tt4.5hi), 
in  R 

xl3  —  Fan  turbine  inlet  slow  response  temperature,  (Tt4.51o), 
in  R 

xl4  —   Fan  turbine  exit  temperature,  (Tt5) ,  in  R 

xl5  —   Duct  exit  temperature,  (Tt6c) ,  in  R 

xl6  —   Duct  exit  temperature,  (Tt7m) ,  in  R 
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pressures  and  11  total  temperatures.   The  inputs  to  the  sys- 
tem are  the  fuel  flow,  the  3  variable  geometry  elements  pre- 
viously discussed,  and  the  customer  bleed  air  percentage. 
The  outputs  include  the  net  thrust,  the  engine  total  airflow, 
turbine  inlet  temperature  and  the  stall  margins.   Addition- 
ally, two  fan  exit  pressure  ratios  are  determined  for  some 
of  the  operating  points.   The  inputs  and  outputs  are  listed 
in  Table  II-2. 

The  technique  used  by  Miller  and  Hackney  perturbed  each 
state,  xi ,  slightly  while  holding  the  other  states  and  inputs 
constant.   This  allows  calculation  of  the  deviations  caused 
by  that  state  in  the  other  states  and  outputs.   These  devia- 
tions then  form  the  F  and  H  matrices  of  equations  (II-5) 
and  (II-6) .   Perturbing  each  input,  ui ,  in  the  steady  state 
configuration  determine  the  G  and  D  matrices.   The  sampling 
time  was  set  at  7  milliseconds. 

The  linear  models  obtained  are  very  good  approximations 
of  the  non-linear  system  but  are  often  too  complex  and/or  do 
not  contain  the  most  convenient  parameterization  to  be  used 
in  a  practical  design.   A  reduction  of  model  order  and,  per- 
haps, augmentation  of  the  linear  model  is  generally  required 
prior  to  using  the  model  for  control  design. 

E.   CONTROL  DEVICE  EFFECTIVENESS 

Before  designing  an  engine  control,  one  must  determine 
the  effect  each  input  has  on  the  system  and  the  responsive- 
ness of  each  input.   The  five  control  inputs  listed  in  Table 
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TABLE  II-2 
Engine  Inputs 
ul   —  Main  burner  fuel  flow,  lb/hr 
u2   —  Nozzle  jet  area,  ft**2 
u3   —   Inlet  guide  vane  position,  deg 
u4   —  High  compressor  variable  vane  position,  deg 
u5   —  Customer  bleed  flow,  % 

Engine  Outputs 

yl  —  Net  thrust,  lb 

y2  —  Total  engine  airflow,  lb/sec 

y3  —  Turbine  inlet  temperature,  R 

y4  —  Fan  stall  margin 

y5  —  Compressor  stall  margin 

y6  —  Fan  exit  delta  p  ratio  (test  data) 

y7  —  Fan  exit  delta  p  ratio  (theory) 
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II-2  are  used  in  controlling  the  F100  engine.   The  customer 
bleed  air  percentage,  u5,  is  essentially  constant  in  steady- 
state  operations  and  can  therefore  be  ignored  in  the  control 
development.   A  study  could  be  made  to  determine  the  merit 
in  extracting  additional  bleed  air  from  the  engine  as  a  pri- 
mary control,  and  dumping  that  portion  not  required  by  the 
customer  systems. 

The  intent  of  this  control  development  is  to  minimize  the 
specific  fuel  consumption  of  the  engine.   At  the  same  time, 
only  small  variations  in  turbine  inlet  temperature,  y3,  and 
total  engine  airflow,  y2 ,  are  desirable  and  no  change  in 
thrust  is  allowable.   One  method  might  employ  using  exhaust 
nozzle  area,  u2 ,  and  fan  inlet  guide  vanes,  u3 ,  inputs  to 
minimize  SFC  at  the  trim  condition  and  then  compensate  for 
thrust  loss  by  using  fuel  to  balance  the  equation. 

Figure  II-5  plots  the  change  in  specific  fuel  consumption 
due  to  incremental  changes  shown  in  Table  II- 3  in  the  inputs 
ul ,  u2 ,  and  u3;  input  u4 ,  the  compressor  variable  guide  vanes, 
exhibit  no  influence  on  SFC.   Proportional  derivatives,  i.e., 
slopes  were  determined  by  making  a  linear  approximation  to 
the  curves.   It  is  clear  in  the  figure  that  nozzle  area  and 
inlet  guide  vane  inputs  (u2  and  u3)  have  the  most  substantial 
effect  on  SFC. 

Figure  II-6  plots  the  net  engine  thrust,  yl,  versus  the 
control  input  changes.   Again  u2  has  the  greatest  effect  and 
ul  and  u3  contribute  about  equally  to  thrust  variations.   It 
is  also  noted  that  a  linear  combination  of  ul  and  u3  can  be 
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TABLE  II -3 
Input  Perturbations 

The  following  perturbations  are  used  throughout  this 
investigation.   They  represent  the  largest  expected  value , 
or  physical  limitation  placed  on  each  control  input. 

Fuel  Flow,  ul  ±2  00  lb/hr 

Nozzle  Area,  u2  ±0.2  sq.  ft. 

Inlet  Guide  Vanes ,  u3  ±4  degrees 

Compressor  Vanes,  u4  ±4  degrees 

Bleed  Air,  u5  0.0  percent 


used  to  offset  the  thrust  change  caused  by  a  change  in  u2 . 
This  is  important  in  meeting  the  zero  net  thrust  change  de- 
sign requirement.  Figure  II-7  shows  the  effect  the  inputs 
have  on  turbine  inlet  temperature,  y3.  Again,  u2  is  the 
dominant  control  input.  These  figures  clearly  demonstrate 
the  need  for  minimizing  the  exhaust  pressure  mismatch  with 
ambient  conditions. 

Table  II-4  is  a  summary  of  the  proportional  derivatives 
extracted  from  the  preceding  figures.   Using  these  derivatives 
the  following  relationships  are  found. 

d(SFC)   -   2.27E-05  d(ul)  +  .0655d(u2) 

-  .00092d(u3)  (II-7) 

d(yl)   =   .9325  d(ul)  -  1351.7d(u2) 

+  18.6275d(u3)  (II-8) 
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TABLE  I I -4 
Summary  of  Derivatives 

Input                    SFC  Thrust  (yl)  FTIT  (y3) 

Fuel  Flow,  ul            2.27E-5  .9325  .12126 

Nozzle  Area,  u2           .0655  -1351.7  283.085 

Inlet  Guide  Vanes,  u3     -.00092  18.6275  -2.1641 

U4  was  found  to  have  no  effect  on  the  parameters  of 
interest . 


d(y3)   =   .12126d(ul)  +  283.085d(u2)  -  2.1641d(u3)    (II-9) 

If  a  decrease  in  u2  of  -.2  square  feet  and  an  increase  of  +4 
degrees  in  u3  are  applied,  a  change  of  -.017  occurs  in  SFC. 
The  fuel  flow,  ul,  is  used  to  balance  the  thrust  equation. 
To  accomplish  this  the  fuel  is  reduced  by  370  pounds  per  hour, 
further  decreasing  specific  fuel  consumption.   The  total  effect 
of  these  three  inputs  at  the  trim  condition  is  a  3.5%  reduc- 
tion in  SFC,  no  thrust  change  and  a  110  degree  decrease  in 
the  turbine  inlet  temperature. 

By  requiring  that  thrust  level  not  change,  the  SFC  per- 
formance is  compared  under  identical  operating  conditions. 
Each  input  can  be  studied  separately  to  determine  its  effect 
on  the  outputs  and  on  SFC.   To  remain  in  the  small  perturba- 
tion regime  inputs  should  not  be  varied  more  than  10%  of  their 
steady-state  values.   Combinations  of  inputs  are  then  chosen 
to  effect  the  most  desirable  improvement  in  SFC  and  minimize 
the  effect  on  engine  performance. 
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This  technique  has  shown  improvement  in  the  specific  fuel 
consumption  at  the  operating  point  chosen  in  the  studies  of 
Reference  1.   The  use  of  the  existing  variable  geometry  to 
improve  the  fuel  efficiency  at  the  high  altitude  cruise 
operating  point  chosen  for  this  investigation  will  be  dis- 
cussed in  a  later  section. 
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III.   CONTROL  DEVELOPMENT 


A.   LINEAR  OPTIMAL  CONTROL  THEORY 
Given  the  linear  system, 


X   =   FX  +  GU  (III-l) 


consider  a  linear  control  law  defined  by 


U  =   -KX  +  UQ  (III-2) 


The  original  system  can  be  modified  to 


X  =   (F  -  GK)X  +  GUQ  (III-3 


The  response  of  the  modified  system,  as  well  as  its  eigen- 
values, is  determined  by  the  effect  of  the  gain  matrix,  K. 
For  a  low  order  system,  the  control  gains  could  be  found  by 
trial-and-error  methods  and  applied  to  the  system  to  determine 
if  the  response  was  acceptable.   This  method  is  impractical 
for  large  order  systems. 

Optimal  control  is  achieved  by  minimizing  a  cost  function, 
J,  which  is  defined  as 

J   =  j  j   [(XTRXX)  +  (UTR2U) ]dt  (III-4) 

where  R,  and  R-  are  symmetric  weighing  matrices;  R,  is  non- 
negative  definite  and  R_  is  positive  definite.   A  third  portion 
can  be  added  to  (III-4)  to  provide  terminal  state  weighing, 
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x  (t1)P1x(t1) 


(III-5) 


where  P,  is  a  non -negative  definite  symmetric  matrix. 

It  is  shown  by  Kwakernaak  and  Sivan  [8],  that  the  optimal 
input,  U.(t) ,  based  on  the  cost  function  is  given  by 


UQ(t)   =   (R2)"1GTP(t) 


(III-6) 


thus  the  knowledge  of  p(t)  solves  the  regulator  problem.   A 
2n  linear  system  is  formed  with  the  optimal  system  behavior, 
x0(t),  and  the  adjoint  variable,  p(t). 


XQ(t) 


P(t) 


-1  T 
-GR2  G1 


,1 


"Rl      "F' 


X  (t) 


p(t) 


(III-7) 


where 


p(t)   -   P(t)XQ(t) 


(III-8 


with 


K(t)   =   R21(t)GT(t)P(t) 


(III-9 


P(t)  is  a  symmetric  non-negative  definite  matrix  that 
satisfies  the  Riccati  equation 
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-P(t)   =   Rx(t)  -  P(t)G(t)R21(t)GT(t)P(t) 

+  P(t)F(t)  +  FT(t)P(t)  (111-10) 


with  the  terminal  condition  of 


P(tL)   =   Px  (III-ll) 


Under  these  conditions,  as  the  terminal  time,  tl,  approaches 
infinity,  the  control  law  asymptotically  approaches  a  unique, 
stable  steady  state  condition.   Thus,  the  matrix  Riccati 
equation,  (111-10) ,  can  be  rewritten 


0   =   R  -PGR~1GTP  +  PF  +  FTP  (111-12) 


and 


K   =   R21gTp  (111-13) 


Many  numerical  methods  are  available  to  solve  the  matrix 
Riccati  equation  for  the  steady-state  control  gains.   Some 
of  these  are  covered  by  Kwakernaak  and  Sivan  [6] .   This 
thesis  employed  a  computer  solution,  0PTSYS4,  developed  by 
Hall  and  Bryson  [7] ,  at  Stanford  University.   A  representative 
output  is  included  in  Appendix  3. 

B.   CONTROL  SYSTEM  MODELS 

At  this  point,  an  analysis  of  control  models  is  made  to 
establish  the  form  of  the  controller  to  be  applied.   A 
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comparison  of  small  and  large  perturbation  controls  must  be 
made.   Take  the  linear  system 

X   =   FX  +  GU  (111-14) 

Y   =   HX  +  DU  (111-15) 

This  system  can  be  made  closed-loop  by  incorporation  of  a 
control  law 

U   =   -  KX  .  (111-16) 

K  can  be  a  time-varying  matrix  or  constant  as  required  by 
the  system.   A  more  general  control  law  would  be 


U   =   U   +  C  (X  -X  )  +  C   A(Y  -Yrt)dt  (111-17) 

m    x     m    '  y      0 


where  A  is  the  diagonal  matrix  whose  diagonal  terms  are  1 
or  0  determining  which  outputs  are  used  in  formulating  the 
input.   C   and  C   are  again  matrices  associated  with  the 
states  and  the  outputs,  respectively.   Additionally  define 


U    =   Un  +  R  t  (111-18) 

m      0    u 


X    =   Xn  +  R  t  (111-19) 

m       Ox 


the  0  subscript  indicates  trim  position/value;  R   and  R   are 
r  F  u      x 

functions  of  u  and  x  associated  with  time.   For  small  pertur- 
bations the  time  varying  terms,  R   and  R  ,  can  be  ignored  and 
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equation  (111-17)  can  be  rewritten 


5U   =   C  5X  +  (C  A6Y  dt  (111-20) 

x     ;  y 


Dropping  the  delta  terminology  and  taking  the  derivative  of 
this  and  using  equation  (111-16) ,  yields 


-  KX   =   C  X  +  C  AY  (111-21) 

x     y 


combining  terms 


(-K-C  )X   =   C  AY  (111-22) 

x        y 


Substituting,  using  the  system  description  of  equations 
(111-14)  and  (111-15)  one  obtains 


(-K-C  ) (FX+GU)   =   C  A(HX+DU)  (111-23) 


and  solving  for  C  , 


C    =   -K-C  A(H-DK) (F-GK)"1  (111-24) 

x         y 


and  using  all  outputs  (A  =  I)  and  equally  weighing  them 
(Cy  =  I)  one  has 


C    =   -K-(H-DK) (F-GK) "1  (111-25) 
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However,  if  the  outputs  are  ignored  (A  =  0),  i.e.,  consider- 
ing the  case  where  small  changes  are  made  in  the  output 
vector,  then 


C   =   -  K  (111-26) 

x 


as  in  the  simplified  control  law  of  equation  (111-16) . 

This  approach  does  not  apply  to  the  large  perturbation 

case.   In  large  perturbation  analysis,  the  time  varying 

terms  cannot  be  ignored.   The  mean  values  of  u  and  x  are  as 

in  equations  (111-18)  and  (111-19) .   R  and  R   are  the  time 
^  u      x 

functions  that  describe  the  path  that  the  mean  values  follow 
over  the  period  of  the  large  perturbation.   The  control  law 
becomes 


5U  -R  t   =   C  (oX  -R  t)  +  /C  A6Y  dt  (111-27) 

u       x      x     ;  y 


Again,    if    all   outputs    are   used  and   equally  weighed,    then 
taking   the  derivative 


-R     -R    t      =      C    X  -C    R      -C    R    t  +C   AY  (111-28) 

uu  xxxxxy 


and   combining    terms 


-C    (F-GK)X+R    C     +R   C    t      =      K(F-GK)X    +    R 
X  X    X  X    X  u 


+    R   t  +  (H-DK) X  (111-29) 
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If  we  set 


R  C   =   R  (111-30) 

x  x      u 


and 


R  C    =   R  (111-31) 

xx      u 


in  effect,  assuming  a  linear  relationship  between  the  states 
and  inputs .   Then 


-C    =   -k  -  (H-DK) (F-GK)"1  (111-32) 

A 


as  in  equation  (111-25)  for  the  small  perturbation  case. 

An  adaptive  control  system  removes  the  need  to  model  the 
large  perturbation  response.   Linear  models  can  be  developed, 
as  described  in  Section  II. E,  to  define  the  large  perturba- 
tion case  as  a  series  of  linear  relations.   The  onboard  digi- 
tal computer  then  schedules  the  control  gains  based  on  the 
linear  models  as  the  flight  conditions  pass  through  the  model 
transition  points.   The  rapidity  of  response  inherent  to  the 
digital  computer  removes  the  requirement  for  modeling  the 
total  dynamic  system,  rather,  only  the  dominant  response  model 
for  the  selected  operating  point.   The  reduced-order  linear 
model  is  far  easier  to  handle  in  formulating  the  control  sys- 
tem and  is  derived  to  match  the  dominant  system  characteristics 
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C.   REDUCTION  OF  LINEAR  MODELS 

DeHoff  and  Hall  [8],  present  a  method  of  reducing  the 
order  of  the  linear  system  based  on  the  dominance  of  the 
states  and  the  crosscoupling  of  the  states  as  observed  through 
eigenvector  analysis.   The  required  order  of  reduction  is 
determined  by  the  number  of  dominant  states  appearing  in  the 
desired  control  bandwidth.   For  the  F100  study,  DeHoff  and 
Hall  chose  a  control  bandwidth  of  0  to  10  rad/sec,  corres- 
ponding to  the  primary  control  device  (fuel  flow  actuator) . 
Their  operating  point  (static,  sea-level,  intermediate  power) 
dictated  that  the  control  be  designed  to  modulate  thrust, 
providing  maximum  thrust  for  take-off  performance.   At  another 
operating  point  a  different  parameter  may  dictate  the  control 
design  and,  thus,  the  primary  control  actuator  and  control 
bandwidth. 

The  method  of  decomposition  involves  finding  the  eigen- 
vector transform  matrix,  T,  such  that 

XT   =   TA  (111-33) 

where  A  is  the  diagonal  eigenvalue  matrix.   Then,  by  defining 
an  alternate  state  vector,  Z,  where 

X   =   TZ  (111-34) 

the  original  system  may  be  written  in  modal  coordinates  as 


Z   =   AZ  +  T~1EU  (111-35) 
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If  a  control  bandwidth  is  determined,  the  eigenvalues 
falling  in  that  range  are  the  only  ones  that  need  to  be  in- 
cluded in  the  initial  reduced  order  model.   The  matrix  may 
be  reordered  by  constructing  a  matrix  Q,  which  has  a  1  in 
the  position  corresponding  to  the  eigenvalue  to  be  included 
(column)  and  the  reordered  position  (row)  of  that  eigenvalue. 
Thus,  if  eigenvalue  5  of  the  original  A  matrix  was  to  be  re- 
ordered to  the  first  position,  a  1  would  be  placed  in  the 
(1,5)  position  of  the  Q  matrix.   Now,  the  reordered  matrix 


A  '   =   QA 


(111-36) 


is  incorporated  in  the  linear  model 


QZ   =   A ' Z  +  QT_1BU 


(111-37) 


letting 


Z  '   =   QZ 


(111-38) 


Z*   =   A  ,Q"1Z'  +  QT_1BU 


(111-39) 


and 


X   =   TQ""1Z* 


(111-40) 


A  matrix  R  can  be  found  such  that 


X»   =   RX   =   TZ ' 


X. 


X. 


(iii-4i; 


42 


This  X'  vector  is  the  reordered  state  vector,  where  the 
elements  contained  in  the  subvector  X..  are  those  states 
associated  with  the  eigenvalues  in  the  control  bandwidth 
and  X2  contains  all  other  states . 

The  linear  system  may  be  rewritten 


Now  define 


and 


X'   =   RFR 


~1X'  + 


RGU 


HR 


"1X'  + 


DU 


(111-42) 
(111-43) 


F1   =   RFR 


-1 


(111-44) 


G'   =   RG 


(111-45) 


H'   =   HR"1 


(111-46) 


11 


21 


12 


22J 


(111-47) 


(111-48) 


H'   -    [H.   i    H0] 
1   l     2. 


(111-49) 
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If  the  states  in  X,  truly  model  the  full  order  system, 
then  X2  will  be  essentially  zero  and 


Xl   =   F11X1  +  F12X2  +  G1U  (111-50) 


0   "   F21X1  +  F22X2  +  G2U  (IH-51) 


Solving  for  X_  in  equation  111-51,  one  has 


X2   '   "  F221f21X1  +  F22lG2U  (III"52) 


This  is  incorporated  into  the  first  equation 


Xl  "   (F11-F12F221f21)X1  +  (GrF12F22lG2)U        (IH-53) 


and  the  output  equation  becomes 


Y   =   (H1-H2F221F21)X1  +  (D-H2F221G2)U  (111-54) 


The  reduced  order  model  can  be  written  as 


where 


X   =   FrX   +  GrU  (111-55) 


Y   =   HrX   +  DrU  (111-56) 


Fr   =   Fli  "  F12F22lF21  (111-57) 
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Gr   =   Gl  -  F12F2llQ2  (HI-58) 


Hr   =   H1-H2F2"1F21  (111-59 


Dr   =   D-H2F221G2  (111-60) 


This  technique  was  used  to  develop  the  program  REDUCR2 , 
documented  in  Appendix  A.   REDUCR2  was  used  to  develop  the 
reduced  order  models  used  in  the  control  development  portion 
of  this  thesis.   Sample  input  and  output  of  REDUCR2  are 
listed  in  Appendix  A. 

D.   REDUCED  ORDER  MODELS 

Three  reduced  order  models  were  generated  with  the  REDUCR2 
computer  program  (Appendix  A) .   The  criteria  for  selecting 
the  desired  states  was  that  used  by  DeHoff  and  Hall,  a  con- 
trol bandwidth  of  0-10  Hertz.   The  eigenvalues  of  the  full 
16th  order  system  were  found  using  the  OPTSYS4  computer  pro- 
gram and  then  associated  to  the  system  states  using  modal 
analysis.   Table  III-l  gives  the  results  of  this  modal  analy- 
sis.  The  five  states  associated  with  the  control  bandwidth 
are  fan  speed  (xl) ,  compressor  speed  (x2) ,  augmentor  pressure 
(x5) ,  burner  exit  temperature  (xlO) ,  and  fan  turbine  inlet 
temperature  (xl3) . 

A  fifth  order  model  was  generated  using  these  states.  It 
was  noted  that  only  three  of  these  states  (the  two  rotational 
speeds  and  the  pressure)  were  dynamically  significant,  i.e., 
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TABLE  III-l 

Operating  Point  Eigenvalues 


State 

Physical  Property 

Eigenvalue 

X4 

P04.5 

-540.23 

X3 

P03 

-156.432 

Xll 

T04 

-49.69  ±  J13.99 

X16 

T07 

X9 

T04Hi 

-  54.12 

X12 

T04.5Hi 

-  42.61 

X7 

T02.5C 

-  22.59 

X14 

T05 

-21.58  ±  J2.51 

XI 5 

T06C 

X6 

T02.5H 

-  18.62 

X8 

T03 

-  14.49 

X5 

P07 

-   7.35 

XI 

Nl 

3.60     control 
band- 

XI 3 

T04.5LO 

-   1.98     width 

X2 

N2 

-   1.31 

X10 

N04LO 

.647 
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TABLE  III -2 
Reduced  Model  Comparisons 


XI 

X2 

X5 

X10 

Xll 

XI 3 

X16 

Yl 

Y3 


Model  #1 

Model  #3 

Model  #2 

Full  Order 

173.51 

175.89 

175.89 

173.41 

125.10 

126.86 

126.86 

125.02 

.326 

.325 

.325 

.326 

4.85 

4.86 

53.51 

53.48 

3.98 

3.98 

23.04 

23.03 

94.48 

174.8 

97.70 

156.99 

11.84 

53.51 

9.3 

53.48 
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represent  energy  storage  in  the  form  of  torque  and  pressure. 
A  third  order  model  was  generated  using  these  three  states 
and  a  comparison  made.   In  this  initial  analysis  these  two 
models  produced  virtually  identical  response  with  respect 
to  the  key  outputs  thrust  (yl)  and  turbine  inlet  temperature 
(y3) ,  see  Table  III-2.   The  0PTSYS4  program  was  used  to  deter- 
mine the  transfer  functions  associated  with  these  two  outputs 
for  each  model.   It  was  found  that  the  fifth  order  model  con- 
tained zeros  near  the  poles  at  -1.9  60  and  -.651,  the  eigen- 
values associated  with  the  two  temperatures  in  the  fifth 
order  model.   This  caused  these  two  poles  to  have  no  effect 
on  the  response  of  the  model.   Thus  the  third  order  system 
sufficiently  models  the  control  bandwidth. 

In  both  models,  the  turbine  inlet  temperature  (y3)  response 
did  not  closely  model  that  of  the  full  order  system.   There- 
fore, a  third  model  was  formed  using  the  three  states  of 
the  third  order  model  plus  the  conjugate  pair  of  poles  repre- 
senting the  turbine  inlet  temperature  (xll)  and  exit  tempera- 
ture (xl6) .   This  model  gave  excellent  thrust  and  temperature 
response  (yl  and  y3 ,  Table  III-2)  due  mainly  to  the  inclusion 
of  the  turbine  inlet  temperature  as  a  system  state.   Table 
III-3  lists  the  three  models,  their  states  and  associated 
eigenvalues . 

Although  the  conjugate  pair  of  eigenvalues  lies  outside 
the  control  bandwidth,  their  inclusion  markedly  improves  the 
output  correspondence  with  respect  to  the  full  order  system, 
see  Table  III-2.   Further  analysis  will  use  the  third  order 
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TABLE  III-3 
Reduced  Model  Eigenvalues 


States 


Model    #1 


XI 


X2 
X5 
X10 


Model    #3 


XI 


X2 


X5 


Xll 


Model    #2 

XI 

X2 

X5 


X13 

XI 6 

-7.175 

-51.65    ± 

-6.947 

Eigen- 

-3.392 

J7.14 

-3.282 

values 

-1.418 

-    1.388 

-1.389 

-  .651 

-    6.909 

-1.960 

-    3.265 
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TABLE  I II -4 

Fifth  Order  Model 


I     t     I    +  + 

OQGQO 

•M  (7s -TOO 

mfNIfM^O'' 
COfNJvOO  co 


OOOO'- 
4-  +  4  4  4 

OQQQQ 

co— «— iOin 

mf-cmo 
ooo-^coo 


oo 
i   i 

OQ 
<^— i 


Oooo 
l    l    I    l 

OQQQ 

cor^OvO 
r^o^co 


OO^Oco-4-cr  vO 


oo  oooo 

4  +  till 

OQ  3QOC 

OO  oooo 

oo  oooo 

0>f  — 'nJ-OcO 

o— «  r— oocr 


I    I    t 

t    1    t 

hMHO-t-4lftH 
1            III 

it       i       ii 

oon^-^ 

— 1  — «<NJi-  <■■  I 

oioo-nr><rin 

cm 

r«-m->tir\ 

ooooo 

ooooo 

ooooooo 

oo 

OOOO 

+•  4>  t  ♦♦ 

♦  ♦»♦<* 

+  i 

♦■  1    I    1    l 

+■  i 

1    1    1    1 

QQGOQ 

ooooo 

oooaooo 

oo 

oooo 

f\lcOO«Of\( 

cor*-ocg>*- 

m-^oir\soorj 

oo 

oooo 

mr*-ir\(\jr» 

^oommr*- 

oorvj  O  o  O  co  <nj 

CO 

oooo 

— *<?  rrir.  co 

<TIT\  vO^^-O 

00(NJOCCOO<f 

— «r\i 

minmrg 

0*00  0000 

ma^coso 

Of~-ocTvomvj- 

0^«— < 

COvfOvj" 

.0<*J— *r^Q- 

mOoo-Oao 

(NvOOm— «t*—  cjn 

0s  T  O  ^u  r«-  ^t-  -O 

I    I 

—      Ill 

</i      l        l    I    1    i 

m-f^ircj^om 
1 

«m— 'Omm 

XHOi   <■   <»■* 

c\jm 

(\jm<Njrvj 

m— i 

(>g^».^rvj 

OOOOO 

—ooooo 

xoo 

oooo 

oo 

OOOO 

♦  4-  4.  4  + 

0C4  +I+  + 

—  4-    I 

i   I    i    i 

4    1 

1    1     1    1 

T^ono 

►-ooooo 

ncaa 

oooo 

oo 

oooo 

tn-o  r-i-o 

<ir\en>4-— <r- 

t—  ^4-cO 

(MOOC0 

C^tn 

NOmiTivT 

r-cr  or»-f- 

3"— i-orfnnif\ 

<mr\l 

vOOrOfNj 

rg<M 

-OT5aom 

CMT\r~-m^<j 

OO—oom 

yi<f\<t 

vO^-oa1 

VIOO 

— »iAO-4" 

:OO^jr\)*l- 

Zmm-*m.J* 

r~in 

rM-^vTi  sj- 

•— 'CJ^O 

mcT1  -^m 

m— 4-f  -o-^ 

OrriajiNO^o 

3  co— « O  »-<  O  r--  o 

o  j-ococ7»-rvO 

•    •     •    •     • 

—•••••• 

O   •  . 

•   •    •   •    • 

X    •    • 

•   •   •   •   • 

rovir*--*— * 

i-ar-— 4-jro 

»-«  — tC*  O  <^  *4f0  — * 

—  -4-  m  o  — •  "O  — « -4" 

I        l    I 

3  1                1 

»—      1 

1         1    1 

QC  1 

in 

CO 

D 

►— 

•—OOftOO 

— .oj  — j^rg-H 

CO—H^- 

mvf  irir^- 

<-4— . 

mmmm 

(^><^>(^\<->  o 

CtOOOOO 

— oo 

OOO1^ 

5TOO 

oooo 

x  ♦  +  I  4  ♦ 

t—  4.  4>  4  4  4- 

or  1    1 

1    1     1     1 

4-    I 

till 

-OQCGQ 

t^OOOOC 

i-ao 

OOOO 

ooo 

oooo 

orr — '<i-t  -o 

•—•■t  mir\  mo 

OOOO^J 

mcor*-f~- 

OiOxO 

vj-mm  — 

>—  IT.— '— *<+\<J 

OO— ivomo 

-^a^rg 

f\4»Avt'vf 

<nTO 

oooo  CO 

<<\j^-m'./*»  — < 

cor\j4-mo 

Ostvf 

— <0"LfMf\ 

3rOr- 

sOCOvOf^ 

7Tr>-'Nmcom 

-JO<M^O  OP* 

h«t 

in^mtM 

SC<ff^ 

— *&^o<t 

»*  ^o  «r  m  ^j 

LJvrr^r--_^rvj 

»—  ocr  o^  <r  o  — • 

U^MOOflD  f\|— < 

►—••••• 

x    •    •    •    •    • 

3    •   • 

•     •     •     •     • 

U_     •    • 

•      •      •      •      • 

^■fM-^rjU^rg 

^-^(^lmm^ 

a_-*c  o— »-«irico 

Qmrgosin-O-^ 

<        III 

2!  1     1    | 

V-        1 

1    1 

UJ    I      I 

1 

-J 

O 

Z> 

UJ 

ao-imoc 

(_JfM-4fVJ— «0 

n  — *fO 

IT  lTMr\  ^0 

KL-Kf 

irv  m  -4*  u"n 

ooooo 

ooooo 

oo 

oooo 

oo 

oooo 

Q  *    1      1*4- 

O   1      1     1*4 

O  1    I 

1   1   1   1 

a  i   i 

till 

UJCQQCQ 

UJQCCQO 

-UOO 

oooo 

LUOO 

oooo 

o^cor-mco 

ocm^"«oir\ 

<_)>*■  >o 

-*<r-oo 

UOff 

— 1-4-— <s0 

ID  J^^jirivrcr 

^sTnT^—CT 

3— <^o 

mor^tn 

^jvoao 

^00  nTvO 

QO^cd-  v0 

Qvj-a^'Moo 

oooo 

mr~cr-^ 

Or-oj 

op--c^in 

UJtr^l^.im^ 

LLU-fOlfiCOO 

UJ<N|sO 

Lnooc*\^- 

ujm<r 

^JOf>-<M 

cr  o^jr*-— «o 

aim  .f-f-  O^cc 

a:CTvrofNicoof\J 

xooo-omoc 

T  '  i 


jjm  }■—•—•—• 
I  I    i  i 


OJ-*v0O-'-J"-iN 

r  i 


OJ  — •  — 'O  T5— •  — «<%J 
Til  I 


50 


TABLE    III-5 
Third  Order  Model 
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model  (model  #2)  and  the  fifth  order  augmented  model  (model 
#3) .   These  models  are  listed  in  Tables  III-4  and  III-5. 

E.   VERIFICATION  OF  REDUCED  ORDER  MODELS 

The  reduced  order  model  must  be  sufficiently  representa- 
tive of  the  original  system  for  the  designed  control  to  pro- 
vide the  desired  full  system  response.   Weinberg  and  Adams 
[9] ,  used  the  F100  non-linear  simulation  to  develop  a  linear 
17th  order  system  at  the  static,  sea-level  point  at  partial 
power  settings.   Their  reductions  to  5th  and  3rd  order  and 
similar  reduction  to  5th  order  by  DeHoff  [10] ,  from  Miller 
and  Hackney's  16th  order  model  show  excellent  correlation  to 
the  non-linear  simulation. 

The  operating  point  chosen  for  consideration  in  this 
thesis  is  30,000  feet,  Mach  =  0.9  and  power  level  angle  of 
67  degrees.   This  corresponds  to  a  near  optimum  cruise  point 
and  is  a  condition  frequently  encountered  in  routine  opera- 
tions.  Section  III.D  lists  the  states  and  associated  eigen- 
values at  this  operating  point  and  discussed  several  reduced 
order  models.   The  dominant  eigenvalues  were  considered  to 
be  those  which  affected  the  models  fast  response  and  had 
dynamic  significance  to  the  engine's  physical  characteristics 
Again,  this  implies  a  control  bandwidth  of  0-10  rad/sec,  how- 
ever, other  models,  using  eigenvalues  outside  the  control 
bandwidth,  were  considered  in  determining  the  best  reduced 
model  for  deriving  the  control  laws. 

One  method  to  analyze  model  correspondence  is  to  use  the 
Bode  plot  across  the  control  bandwidth.   Figure  III-l  is  a 
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Bode  comparison  of  the  DeHoff  reduced  model  and  the  Miller 
and  Hackney  16th  order  model  at  static,  sea-level,  inter- 
mediate power  for  the  0  to  10  rad/sec  control  bandwidth. 
It  is  clear  that  through  the  0-10  Hertz  bandwidth  the  two 
models  show  very  close  correspondence;  this  similarity  breaks 
down  in  the  higher  frequency  region,  but  the  response  time  of 
the  eigenvalues  outside  the  control  bandwidth  is  considered 
to  be  too  short  to  affect  transient  behavior.   The  phase 
response,  not  shown,  follows  the  same  trends. 
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IV.   ANALYSIS  OF  CONTROL  SYSTEMS 

A.   COST  FUNCTION  WEIGHING  FACTORS 

The  control  law  for  the  engine  regulating  system, 


U   =   Um  +  Cx(X_Xm)  +  /  CyA(Y"V  dt  (IV"X) 


can  be  reduced  as  shown  in  Section  III.B  to 


U   =   U   -  KX  (IV-2) 

m 


when  considering  small  perturbations  and  ignoring  the  out- 
puts from  the  system.   In  most  circumstances,  the  U  ,  the 
trim  condition,  would  be  considered  as  zero  for  steady  state 
analysis.   The  linear  models  developed  by  Miller  and  Hackney 
[5],  set  all  inputs  to  zero  at  the  steady-state  operating  point 
for  which  the  model  is  defined. 

The  K  matrix  is  the  solution  to  the  optimum  regulator 
problem  discussed  in  Section  III. A.  The  0PTSYS4  computer 
solution  used  in  this  thesis  uses  the  cost  function 


J   =  \  I   [YTA  Y  +  UTR0U]dt  (IV- 3 

2  j  y       2 


where  A   is  the  weighting  matrix  for  the  output  vector  and 
replaces  the  need  for  the  R,  matrix  in  equation  (III-7)  . 

Kirk  [11],  discusses  the  effect  that  weighting  matrices 
have  on  system  performance.   In  general,  the  larger  the 
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magnitude  of  the  matrix  element  the  faster  its  associated 
parameter  (system  state,  input  or  output)  stabilizes.   An 
additional  factor  in  determining  the  magnitude  of  the  diagonal 
matrix  elements  is  the  magnitude  of  their  associated  param- 
eters.  Given  the  cost  function 


1  t      T 
J   =  j  j    X  RjX  dt 


(IV-4) 


where 


R. 


1_ 
2 


l_ 
2 


X   = 


xl 


x2 


(IV-5) 


The  cost  function  then  becomes 


j  •  |    [(^)2  +  (^-)2]  dt 


(IV-6) 


The  desire  is  to  minimize  the  cost  function,  J,  with 
respect  to  time.   Thus,  taking  the  time  derivative  and 
setting  it  equal  to  zero 


a       b 


=   0 


(IV-7) 


This  is  the  solution  which  produces  a  minimum  cost.   The 
values  of  a  and  b  are  the  largest  expected  values  of  xl  and 
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x2 ,  respectively,  and  represent  an  equal  effect  on  the  cost 
function  by  each  state.   Should  either  a  or  b  be  decreased, 
the  optimal  control  developed  would  have  to  increase  the 
control  input  to  the  other  state  to  minimize  the  cost  func- 
tion.  A  nominal  fluctuation  limit  of  ±5%  will  be  considered 
as  a  small  perturbation  in  the  initial  control  analysis. 

Designing  a  control  to  minimize  the  thrust  specific  fuel 
consumption  (SFC)  of  the  engine  requires  limiting  the  thrust 
and  the  fuel  flow  fluctuations.   The  weighting  matrix  ele- 
ments corresponding  to  these  two  parameters  are,  therefore, 
to  be  increased  with  respect  to  the  other  parameters.   Looking 
at  only  the  small  perturbation  response,  it  is  noted  that 
very  little  change  in  the  total  engine  airflow  (y3)  is  ex- 
pected and,  as  before,  there  is  no  change  in  customer  bleed 
air  requirements  (u5) .   These  two  parameters  can  be  neglected, 
as  their  response  characteristics  have  no  effect  on  the  desired 
SFC  response.   To  demonstrate  the  effect  of  the  different  ele- 
ment magnitudes  on  the  control  gain  matrix,  the  0PTSYS4  pro- 
gram was  run  for  the  3  cases  delineated  in  Table  IV-1.   Case 
I  was  developed  to  allow  each  output  and  input,  with  the  excep- 
tion of  the  two  pressure  ratios  (y6  and  y7)  and  the  bleed  air 
percentage  (u5) ,  to  effect  equal  contribution  to  the  cost 

function.   The  diagonal  elements  of  both  A   and  R-  normalize 

y      2 

the  squared  value  to  the  order  of  one.  Case  II  uses  all  of 
the  same  elements,  except  the  thrust  (yl)  and  the  fuel  flow 
(ul)  are  weighted  more  heavily,  to  the  magnitude  of  ten.   This 
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TABLE  IV-1 
Case  Definition 

A  and  R?  are  diagonal  matrices,  only  the  significant 
elements  are  listed. 


Case  I 


A   =   (.0001   1   .0001   100   100   le-8   le-8) 


R2   =   (.0001   100   .01    .01  le-8) 


Case  II 
A   =   (.01   1   .0001   100   100   le-8   le-8) 
R2   =   (.01   100   .01    .01  le-8) 


Case  III 

A   =   (1   1   .0001   100   100   le-8   le-8) 

Y 

R2   =   (1   100   .01    .01  le-8) 


heavier  weighting  should  reduce  the  fuel  input  to  the  system 
and  moderate  the  thrust  fluctuations  as  the  control  returns 
the  system  to  the  steady  state  condition.   Case  III  increases 
the  weighting  of  yl  and  ul  to  achieve  a  normalized  magnitude 
of  100.   This  extremely  heavy  weighting  should  eliminate  fuel 
from  the  control  and  reduce  thrust  variations  to  absolute 
minimum.   The  whole  intention  of  the  controller  design  is  to 
damp  out  state  perturbations  without  using  fuel  or  altering 
the  thrust  level,  i.e.,  making  no  change  in  specific  fuel 
consumption  to  effect  control. 
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These  three  cases  were  implemented  on  the  two  reduced 
order  models  and  simulations  run  for  a  perturbation  of  -50 
RPM  in  fan  speed  (xl) .   The  fifth  order  model  returned  the 
fan  speed  to  95%  of  steady  state  in  .5  seconds  using  Case  I, 
but  increased  fuel  flow  by  10  pounds  per  hour  initially  and 
required  an  increase  in  fuel  consumption  of  .2%  over  the 
half  second  period.   For  Case  II,  the  fifth  order  model  re- 
turned to  9  5%  of  the  steady  state  fan  speed  in  approximately 
.3  seconds,  while  using  almost  8  0%  less  fuel  than  the  Case  I 
controller.   This  8  0%  reduction  in  fuel  consumption  indicates 
that  the  increased  weighting  of  Case  II  has  accomplished 
exactly  what  was  desired.   Case  III  was  run  for  the  fifth 
order  model  and,  although  virtually  no  fuel  was  used  in  the 
control,  the  system  experienced  oscillations  in  fan  speed  of 
±8  RPM  through  the  3  second  simulation  run.   Additionally, 
the  Case  III  controller  caused  rapid  changes  of  ±4  degrees  in 
the  compressor  guide  vanes  (u4) ,  the  regulatory  limit  imposed 
for  the  simulation,  and  also  large  fluctuations  in  the  inlet 
guide  vanes  (u3) .   One  would  achieve  results  of  similar  magni- 
tude for  each  of  the  cases  assuming  a  +1%  perturbation  of  all 
states  in  the  reduced-order  model.   Case  I  is  too  fuel  de- 
pendent and  causes  large  changes  in  u3 ,  Case  III  also  requires 
large  changes  in  U3  that  could  cause  a  slew  rate  request  larger 
than  that  available  from  the  control  actuator. 

In  modeling  the  control,  the  actuator  and  physical  limi- 
tations of  each  control  must  be  considered  and  appropriate 
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limitations  placed  on  the  control  law  to  avoid  overshoots. 
Table  IV-2  gives  the  actuator  limits  and  maximum  rates 
applicable  to  the  F100  engine  [4].   A  stepping  limit  of  .01 
seconds  is  assumed  in  determining  the  rate  limits  and  fuel 
flow  fluctuations  are  limited  to  ±10%  of  the  steady-state 
value.   No  compensation  is  allowed  for  hysteresis  or  non- 
linear dynamics  of  the  actuators;  these  are  considered 
minimal  in  the  small  perturbation  case,  but  would  need  to 
be  incorporated  for  large  perturbation  control  modeling. 

TABLE  IV-2 
Actuator  Limitations 

Input  Maximum  Minimum  Rate  Limit 

Fuel  flow  16300  lb/hr  450  lb/hr  15800  lb/hr/sec 

Nozzle  area  6.4  sq.  ft.    2.8  sq.  ft.  3.6  sq.ft. /sec 

Inlet  guide  vanes      0  deg.         -40  deg.       4  8  deg/sec 
Compressor  vanes       4  deg.         -40  deg.       4  0  deg/sec 

It  is  evident  that  the  Case  III  weighting  is  not  accepta- 
ble because  of  the  oscillations  and  the  rapid,  large  amplitude 
fluctuations  in  the  guide  vane  positions  (u3  and  u4) .   Case  I 
is  also  inappropriate  for  the  specific  fuel  consumption  mini- 
mizing control  since  it  uses  fuel  as  the  primary  control  in- 
put.  For  these  reasons,  Case  II  is  chosen  as  the  weighting 
for  all  control  modeling  to  be  implemented  in  later  sections. 

Once  the  integral  portion  of  equation  (IV-1)  and  the  feed- 
back gains  are  determined,  the  constant  portion,  U  ,  must  be 
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analyzed.   In  most  situations  this  term  would  be  neglected. 

The  steady  state  analysis  is  now  performed  to  determine  if 

the  specific  fuel  consumption  can  be  reduced  by  making  small 

changes  to  the  trim  position  of  the  variable  geometry.   An 

analysis  of  control  changes  versus  SFC,   as  was  done  at  the 

static  sea-level  condition  in  Section  II. D,  must  again  be 

conducted.   A  minimum  SFC  obtainable  relates  to  certain  input 

values  in  U  which  is  incorporated  into  the  control  law,  im- 
m  r 

proving  the  steady-state  performance  and  controlling  the 
fluctuations  of  the  system  which  will  have  the  greatest  effect 
on  the  SFC. 

As  was  done  at  the  static  sea-level  operating  point,  a 
small  perturbation  in  each  input  was  made  and  its  effect  on 
the  SFC,  thrust  and  other  engine  conditions  tabulated.  The 
goal  of  these  analyses  is  to  develop  a  trim  input,  U  ,  that 
improves  the  fuel  economy  at  steady  state  without  causing 
significant  change  in  the  operating  conditions  (both  states 
and  outputs) . 

The  effects  of  the  control  inputs  on  several  parameters 
are  plotted  on  Figures  IV-1  through  IV-4 .   Figures  IV-1 
shows,  again,  the  large  effect  nozzle  area  (u2)  and  fuel  flow 
(ul)  have  on  specific  fuel  consumption,  the  effect  of  inlet 
guide  vanes  (u3)  is  drastically  reduced  from  that  observed 
at  static,  sea  level  operation.   Figure  IV-2  indicates  the 
offsetting  effect  of  ul  and  u2  and  again  indicates  a  zero 
net  thrust  control  can  be  found.   Fan  speed  (xl)  and  compressor 
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speed  (x2)  effects  are  shown  in  Figures  IV- 3  and  IV- 4,  respec- 
tively.  At  this  operating  point  the  inlet  guide  vane  posi- 
tion (u3)  may  be  most  useful  in  controlling  the  fan  speed 
fluctuation  and  thereby  stabilizing  the  fan  stability  margin 
(y4) .   A  sample  control  would  close  u2  by  .2  square  feet  and 
reduce  fuel  flow  by  57  pounds  per  hour.   This  combination 
produces  a  1.4%  decrease  in  specific  fuel  consumption  and 
results  in  a  4%  reduction  in  fan  speed.   A  trim  condition 

control  input,  U  ,  would  be 

m 


U 


m 


r-57" 

-.2 

0 

0 

L  o    J 

B.   REDUCED  ORDER  SIMULATION  ANALYSIS 

The  model  developed  in  Section  III.C  and  Section  III.D  for 
the  high  altitude  flight  operating  point  was  used  to  initially 
analyze  the  control  law  developed  in  Section  IV. A.   A  reduced 
order  controller  is  simpler  and  cheaper  to  implement  through 
a  digital  computer.   It  requires  less  computer  space  and  far 
less  computing  time  than  the  full-order  linear  model  developed 
controller. 

A  comparison  of  the  three  cases  for  which  gain  matrices 
were  found  indicated  that  Case  2  gave  the  most  desirable  re- 
sults.  The  Case  2  weighting  matrix  was  applied  to  the  third 
order  reduced  system  and  perturbations  applied  to  various 
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states.   The  same  weighting  matrix  was  applied  to  the  fifth 
order  model  (#3)  and  similar  perturbations  applied.   Figure 
IV-5  shows  the  system  response  to  a  50  RPM  perturbation  of 
the  fan  speed.   The  fifth  order  system  has  a  settling  time 
of  approximately  .45  seconds  whereas  the  third  order  system 
has  one  of  1.5  seconds.   This  variation  is  accounted  for  by 
the  two  fast  response  poles  in  the  fifth  order  model, 
(-51.65, j7. 14) . 

These  results  must  be  compared  to  the  16th  order  system 
to  determine  which  model  most  closely  approximates  the  full 
order  system  response. 

C.   FULL  ORDER  SIMULATION  ANALYSIS 

The  proof  of  the  viability  of  a  reduced  order  controller 
is  its  ability  to  control  the  full  order  system  in  the  same 
manner  as  predicted  in  the  reduced  order  model.   The  two 
reduced  models  analyzed  in  Section  IV. B  were  applied  to  the 
sixteenth  order  linear  model.   Additionally,  the  Case  2  weight- 
ing matrix  was  applied  and  a  full  order  controller  was  found 
and  implemented.   All  simulation  was  done  using  the  IBM  Con- 
tinuous System  Modelling  Program,  Model  III. 

All  three  models  were  compared  with  similar  perturbations 
applied.   Figure  IV-6  compares  the  full  order  and  fifth  order 
response  to  the  50  RPM  speed  perturbation.   Figure  IV-7  com- 
pares the  full  order  and  third  order  response  to  the  same 
perturbation.   Figure  IV-8  compares  the  two  reduced  order 
controllers  and  the  full  order  controller  for  the  same  fan 
speed  perturbation. 
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The  results  show  satisfactory  agreement  between  the  re- 
duced order  and  full  order  implementation  of  the  reduced 
order  controller.   This  proves  the  viability  of  the  con- 
troller and  implies  that  comparable  results  will  be  found 
if  the  control  is  implemented  in  the  non-linear  dynamic 
system.   The  fifth  order  controller,  with  two  poles  outside 
the  control  bandwidth,  most  closely  models  the  response  of 
the  full  order  controller.   As  was  seen  previously,  the 
inclusion  of  the  conjugate  pair  of  poles  dramatically  im- 
proves the  response  time  of  the  model  and  provides  a  much 
closer  approximation  to  the  full  order  system  response. 

This  indicates  that  the  fifth  order  model  is  sufficient 
for  the  control  development  process  at  this  operating  point 
The  two  fast  response  poles  outside  the  control  bandwidth 
must  be  included  to  achieve  this  close  simulation. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  study  of  implementing  a  reduced  order  controller 
that  is  flight  regime  dependent  shows  that  at  a  specific 
operating  point  an  improvement  in  specific  fuel  consumption 
can  be  made.   The  conclusions  reached  at  the  end  of  this 
research  are: 

1.  The  variable  geometry  existing  in  the  state-of-the- 
art  turbofan  engines  can  be  used  to  further  improve  the 
specific  fuel  consumption; 

2.  A  reduced  order  controller  can  be  found  and  imple- 
mented that  concentrates  control  authority  on  the  variable 
geometry  rather  than  the  fuel  control,  thus  furthering  the 
fuel  savings;  and 

3.  A  similar  controller  can  be  found  and  implemented  at 
any  given  operating  point  throughout  the  flight  regime,  these 
controllers  could  then  be  organized  into  a  single  digital 
machine  to  provide  a  regime-dependent  engine  control  to  im- 
prove the  airframe/engine  interaction  in  all  flight  conditions 

B.  RECOMMENDATIONS  FOR  FUTURE  STUDY 

One  obvious  extension  of  this  work  is  to  find  further 
reduced  order  controllers. 

Additionally,  work  can  be  done  to  fully  implement  these 
reduced  order  controllers  in  the  non-linear  model  using 
either  the  Pratt  &  Whitney  transient  simulation  digital  model 
or  the  hybrid  computer  model. 
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The  effect  on  fatigue  life  and  life  cycle  cost  due  to 
the  "fuel  saving"  control  schemes  presented  was  not  deter- 
mined.  It  is  therefore  recommended  that  these  factors  be 
explored  in  order  to  fully  analyze  the  usefulness  of  these 
control  schemes. 
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